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Abstract: Measurements at 100 TeV and above are an important goal for the next generation of high 
energy 7-ray astronomy experiments to solve the still open problem of the origin of galactic cosmic rays. 
The most natural experimental solution to detect very low radiation fluxes is provided by the Extensive 
Air Shower (EAS) arrays. They benefit from a close to 90% duty cycle and a very large field of view 
(—2 sr), but the sensitivity is limited by their angular resolution and their poor cosmic ray background 
discrimination. Above 10 TeV the standard technique for rejecting the hadronic background consists in 
looking for "muon-poor" showers. 

In this paper we discuss the capability of a large muon detector (A;j=2500 m^) operated with an EAS 
array at very high altitude (>4000 m a.s.l.) to detect 7-ray fluxes around 100 TeV. Simulation-based 
estimates of energy ranges and sensitivities are presented. 
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Introduction 

The recent TeV results of the HESS experiment 
suggest the existence of a population of galactic 
7-ray sources whose emission extends beyond 10 
TeV in the 5 to 15% of Crab flux range (for E > 1 
TeV). These sources, associated with nearby shell- 
type or plerionic SNRs, the most probable facto- 
ries of galactic cosmic rays, can be studied detect- 
ing gamma-rays (and neutrinos) emission in the 
VHE/UHE energy domain. Therefore, a detector 
capable to perform a continuous all-sky survey at a 
level of about a percent of the Crab flux up to 100 
TeV is needed. The search and study of "Cosmic 
PeVatrons " and their surrounding regions is one of 
the main scientific issues to be addressed by the 
next generation of ground-based 7-ray astronomy 
detectors [1]. 

Current experiments are not able to reach 100 TeV 
because their limited collection area makes the re- 
quired exposures too long. Extrapolating the galac- 
tic source spectra measured by HESS, it appears 
that future experiments need to achieve at least 
100 km^ h in order to obtain meaningful measure- 
ments at 100 TeV. The most natural experimen- 



tal solution that provides such a large exposure is 
given by the EAS arrays observing each source for 
— 1500 h/year As an example, an experiment like 
ARGO-YBJ already results in an exposure of —15 
km^ -h/year, with an angular resolution (—0.2° at 
10 TeV) near the best value attainable by a sam- 
pling array. In addition, the EAS arrays are the 
only ground-based detectors allowing simultane- 
ous and continuous coverage of a significant frac- 
tion of the sky (about all that overhead). Their 
large field of view and high duty cycle (>90 %) 
suit to perform a 7-ray sources population sur- 
vey at VHE/UHE energies. But more important 
is their unique potential that allows to have an ef- 
fective monitoring of the 7-ray activity of a large 
number of highly variable sources like blazars and 
microquasars, as well as the possibility of inde- 
pendent detection and study of GRBs [1]. In ad- 
dition, the recent observations of unidentified ex- 
tended sources from the Galactic plane [2] and in 
the Cygnus region [3] reported by the Milagro Col- 
laboration demonstrate the strength of EAS arrays 
in finding diffuse and extended sources. Therefore, 
the discovery science could be a feature of EAS ar- 
rays. 
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The limited sensitivity in detecting 7-ray point 
sources, characteristic of EAS experiments, is 
mainly due to their poor gamma/hadron separa- 
tion power, limited angular resolution and high 
energy threshold. Exploiting a full coverage ap- 
proach at very high altitude leads to the improve- 
ment of the angular resolution to ^^0.2° and to the 
reduction of the energy threshold well below the 
TeV region. The standard technique to perform a 
gamma/hadron discrimination above 10 TeV with 
EAS arrays consists in looking for "muon-poor" 
showers. 

In this paper we discuss the capability of a large 
muon detector (A^=2500 m^) operated with an 
EAS array at very high altitude (>4000 m a.s.l.) 
to detect 7-ray fluxes up to 100 TeV. An estima- 
tion based on the simulation of energy ranges and 
sensitivities is reported. 

7-hadron separation 

We have simulated, via the Corsika/QGSJet code 
[4], proton- and 7-induced events with energy 
spectra [5, 6] ranging from 0. 1 TeV to 1 PeV at an 
observational level of 4300 m asl (606 g/cm^), cor- 
responding to the YangBaJing Cosmic Ray Labo- 
ratory (Tibet, P.R. China). Two EAS arrays there 
are located: the conventional Tibet AS7 experi- 
ment and the full coverage detector ARGO-YBJ. 
From the study of the shower phenomenology it 
results: 

( 1 ) The muon content of proton-induced show- 
ers exceeds, as expected, that of 7 showers 
by about two orders of magnitude. 

(2) The muon component in 7-induced events 
has a flatter lateral distribution for core dis- 
tances smaller than a few tens of meters. 

(3) The number of low energy e.m. particles 
(e^, 7) in 7 showers exceeds the correspond- 
ing number in proton showers (see Fig. 1). 

(4) Inside the core region (R<20 m), for particle 
energies above about 400 MeV, the number 
of e='= and 7 in proton showers is higher that 
in 7 showers. 

(5) Outside the core region (R>20 m) the num- 
ber of e.m. particles in proton showers starts 
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Figure 1: Integral energy spectra of electromag- 
netic particles (e='=, 7) produced at 4300 m asl by 
proton- (dotted red lines) and 7-induced (full black 
lines) showers in two different distance ranges 
from the core. 

exceeding that in 7 showers at lower ener- 
gies (about 100 MeV). For E>1 GeV the 
two components differ by more than 1 or- 
der of magnitude: high energy e.m. particles 
outside the core strongly indicate proton- 
induced showers. This number is compara- 
ble to the number of muons. 

In addition, the high energy tail of secondary e.m. 
particles in proton-induced showers is responsi- 
ble for non-uniformity in the spatial distribution 
of shower particles. Consequently, an appropriate 
muon detector operated outside the shower core re- 
gion (R>30 m), evaluating all available high en- 
ergy (E> 300 MeV) secondary components of air 
showers {pL, e^, 7), allows a 7/hadron discrimi- 
nation which goes considerably beyond ordinary 
muon counting [7]. The layout for the muon de- 
tector investigated in this preliminary study con- 
sists of 4 detectors each large 24x26 m^ (total area 
2500 m^) simmetrically distributed on the bound- 
ary of a 150 X 150 m^ large EAS array (for exam- 
ple, out of the four sides of the main building of the 
ARGO-YBJ experiment [8]). The envisaged de- 
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Figure 2: Q-f actor vs. minimum multiplicity for 
different cuts (2, 3, 4, 5 /i). Only muon detectors 
outside the shower core region (R> 30 m) have 
been taken into account. 



tector is constituted by 5 tracking planes shielded 
by ^1 m of concrete in order to absorb particles 
with energy lower than w400 MeV. Therefore, in 
the following we will refer to "muons" as to all the 
high energy particles. 

Compared to hadronic showers, the muon content 
distribution of 7-induced events is significantly 
shifted toward lower values of and includes a 
much larger fraction of events with no muons de- 
tected. Therefore, we define as muon-poor those 
events whose is less than a value chosen to op- 
timize the sensitivity to 7-ray sources. In order to 
evaluate the improvement in sensitivity due to the 
presence of a muon detector, we calculate the so- 
called "Q-factor value" Qf — e^(rp)/^l — Vp, 
where Vp is the fraction of cosmic ray background 
rejected with a given cut and e~f{rp) is the frac- 
tion of 7 showers which survives after this cut. 

Two different analysis procedures have been per- 
formed: (1) all muon detectors have been used 
in the calculation of the muon content; (2) only 
muon detectors 30 m farther from the recon- 
structed shower core position have been taken into 
account. In Fig. 2 the Q-factor values as a func- 



tion of multiplicity for different cuts are shown. 
Only /i-detectors outside the shower core region 
(R> 30 m) have been taken into account to cal- 
culate the muon content. The multiplicity is the 
number of charged particles Nch sampled by a 
lOOx 100 m^ full coverage detector whose 4 sides 
are surrounded by the 4 /i-detectors. 

Sensitivity to 7-ray sources 

In order to evaluate the Minimum Detectable Flux 
(MDF) improvement due to the calculated Q- 
factors, we considered as a reference sensitivity the 
ARGO-YBJ MDF calculated for internal events, at 
5 standard deviations in one year of observation for 
a Crab-like source [8] (red curve of Fig. 3). The 
results of these calculations with procedure (2) are 
shown in Fig. 3 by the blue curve. For a conserva- 
tive estimation of the MDF we exclude cuts at very 
low (0, 1 and 2) muon content, where there is a po- 
tential background from mismatched events, and 
we consider that only the showers with <3 muons 
on the detectors are due to a 7-ray. At 30 (80) TeV 
the improvement in the MDF, due to the rejection 
of showers with > 3 detected muons, is a factor 25 
(100), because at sufficiently high energy the 7-ray 
measurement is background-free. 

The rejection power for dilTuse 7-rays 

The detection of an isotropic diffuse flux of UHE 
photons depends not only on the size and the so- 
phistication of the detectors but also on the proper- 
ties of the 7-ray showers. In fact, if these events 
have the same muon content as hadronic show- 
ers the detection of any diffuse 7-ray fluxes would 
be impossible. In order to evaluate the rejection 
power for proton-induced showers we have studied 
how frequently hadronic showers fluctuate in such 
a way to have a low muon content indistinguish- 
able from 7-induced events. In Fig. 4 we show the 
probability that a photon or proton shower gener- 
ates a given number of muons on the simulated de- 
tector, assuming for both primaries the same spec- 
tral index -2.7. We note that the fluctuations about 
the mean values do not follow the Poisson statis- 
tics. The inclusive probability distributions are re- 
ported in the lower panel of the figure. For pro- 
ton (photon) showers the histograms bars in the in- 



Gamma-Ray Astronomy around 100 TeV 




E(Ge V) Number of muons above 1 GeV 



Figure 3: Minimum Detectable Flux (1 year, 5a) 
compared with the Crab spectrum (black line). The 
red curve shows the ARGO-YBJ sensitivity, the 
blue one the MDF due to a 2500 m^ /i-detector 



Figure 4: (a) Probability that a photon or proton 
shower generates a given number of muons. (b) In- 
clusive probabiUty for the histogram shown in (a). 



<3 <4 <5 

Fraction of 7-ray 

retained 85% 94% 97% 

Background level IQ-'^ 6-10^5 
Table 1: 

elusive distributions represent the probabiUty that a 
shower contains less (more) muons than the upper 
(lower) limit of the bin. The plot refers to showers 
with Nch > 13000 (photon median energy ~150 
TeV). We have conservatively fitted the proton dis- 
tribution overestimating the background in order to 
predict their level at very low muon content [9]. 
In Table 1 we list the implications of 7-ray detec- 
tion for showers with Nch > 13000. In this en- 
ergy region a significant fraction of the cosmic rays 
could be heavy nuclei, yielding larger muon num- 
bers compared to proton showers alone. There- 
fore, we are studying here a worst-case scenario, 
as protons are most likely to fake 7-ray s. From 
MC calculations we conclude that diffuse 7-rays 
of 150 TeV energy can be observed down to a 
level ^10~^ of the cosmic ray background if the 
systematic uncertainties of the detector are under- 
stood at the same level. 



Conclusions 

Recent discoveries in the TeV 7-ray astronomy 
provide strong motivations for extending the mea- 
surements to energies >100 TeV where 7-ray 
emission models can be better discriminated. HAS 
arrays naturally provide large exposures allowing 
simultaneous and continuous coverage of ^2 sr 
without facing serious technological challenges. A 
suitable large muon detector operated with an HAS 
array at very high altitude permits a 7 measurement 
essentially background-free for E >50 TeV. 
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